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I. INTRODUCTION 

Measurements of time-dependent CP asymmetries in 
B° meson decays through b — > ccs amplitudes have pro- 
vided crucial tests of the mechanism of CP violation in 
the Standard Model (SM) [1]. These amplitudes con- 
tain the leading 6-quark couplings, given by the Cabibbo- 
Kobayashi-Maskawa [2] (CKM) flavor mixing matrix, for 
kinematically allowed transitions. Decays to charmless 
final states such as (f)K , tt K , rj'K , ljK°, K + K~K°, 
/o(980).ff are CKM-suppressed b — > qqs (q = u,d,s) 
processes dominated by a single loop (penguin) ampli- 
tude. This amplitude has the same weak phase [3 = 
arg (-V c dV* b /VtdV t * b ) of the CKM mixing matrix as that 
measured in the b — > ccs transition, but is sensitive to the 
possible presence of new heavy particles in the loop [3]. 
Due to the different non-perturbative strong-interaction 
properties of the various penguin decays, the effect of 
new physics is expected to be channel dependent. 

The CKM phase (3 is accessible experimentally through 
interference between the direct decay of the B meson 
to a CP eigenstate and B°B° mixing followed by decay 
to the same final state. This interference is observable 
through the time evolution of the decay. In the present 
study, we reconstruct one B° from T(4S) — ► B a B°, which 
decays to the CP eigenstate wK®, rj'K , rfK^, or tt°K° 
(Bcp)- From the remaining particles in the event we also 
reconstruct the decay vertex of the other B meson (Stag) 
and identify its flavor. The distribution of the difference 
At = tcp — ttag °f the proper decay times tcp and itag 
of these mesons is given by 

e -|At|/r 

/(A*) = ~ ^{1± (1) 
[-rjfSf sin(Am d Ai) - C f cos(Am d At)]} 

where rjf is the CP eigenvalue of final state / (—1 for 
luK° s , r,'K° s , and tt°K° s ; +1 for rj'K®). The upper (lower) 
sign denotes a decay accompanied by a B° (B°) tag, r is 
the mean B° lifetime, and Ara^ is the mixing frequency. 

A nonzero value of the parameter Cf would indicate 
direct CP violation. In these modes we expect Cf = 
and —rjfSf — sin2/3, assuming penguin dominance of the 
b — > s transition and neglecting other CKM-suppressed 
amplitudes with a different weak phase. However, these 
CKM-suppressed amplitudes and the color-suppressed 
tree diagram introduce additional weak phases whose 
contributions may not be negligible [4-7]. As a con- 
sequence, the measured Sf may differ from sin2/3 even 
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within the SM. This deviation ASf = Sf — sin2/3 is esti- 
mated in several theoretical approaches: QCD factoriza- 
tion (QCDF) [4, 8], QCDF with modeled rescattering [9], 
soft collinear effective theory (SCET) [10], and SU(3) 
symmetry [5, 7, 11]. The estimates are channel depen- 
dent. Estimates of AS from QCDF are in the ranges 
(0.0,0.2), (-0.03,0.03), and (0.01,0.12) for u>K°, rj'K , 
and ir°K®, respectively [8, 10, 12]; SU(3) symmetry pro- 
vides bounds of (-0.05,0.09) for rj'K and (-0.06,0.12) 
for 7r°i\T° [11]. Predictions that use isospin symmetry to 
relate several amplitudes, including the J = § B — > Kn 
amplitude, give an expected value for S n o K o near 1.0 in- 
stead of sin2/3 [13]. 

We present updated measurements of mixing-induced 
CP violation in the B° decay modes uK°, rj'K , and 
n K°, which supersede our previous measurements [14- 
16]. Significant changes to previous analyses include 
twice as much data for ojK®, 20% more data for rj'K 
and ir°K°, improved track reconstruction, and an addi- 
tional decay channel in rj'K . Despite the modest in- 
crease in data, the uncertainties on S v i K o and C v i K o de- 
crease by 20% and 25%, respectively. Measurements in 
these modes have also been made by the Belle Collabo- 
ration [17, 18]. 



II. THE BABAR DETECTOR AND DATASET 

The results presented in this paper are based on 
data collected with the BABAR detector at the PEP-II 
asymmetric-energy e+e~ storage ring, operating at the 
Stanford Linear Accelerator Center. At PEP-II, 9.0 GeV 
electrons collide with 3.1 GeV positrons to yield a center- 
of-mass energy of ^s = 10.58 GeV, which corresponds to 
the mass of the T(4S) resonance. The asymmetric en- 
ergies result in a boost from the laboratory to the e + e~ 
center-of-mass (CM) frame of /?7 ss 0.56. We analyze the 
entire BABAR dataset collected at the T(4S') resonance, 
corresponding to an integrated luminosity of 426 fb _1 
and (467 ± 5) x 10 6 BB pairs. We use an additional 
44 fb _1 of data recorded about 40 MeV below this energy 
(off-peak) for the study of the non-B°B° background. 

A detailed description of the BABAR detector can be 
found elsewhere [19]. Surrounding the interaction point 
is a five-layer double-sided silicon vertex tracker (SVT) 
that provides precision measurements near the collision 
point of charged particle tracks in the planes transverse 
to and along the beam direction. A 40-layer drift cham- 
ber (DCH) surrounds the SVT. Both of these tracking 
devices operate in the 1.5 T magnetic field of a super- 
conducting solenoid to provide measurements of the mo- 
menta of charged particles. Charged hadron identifica- 
tion is achieved through measurements of particle energy 
loss in the tracking system and the Chcrcnkov angle ob- 
tained from a detector of internally reflected Chcrcnkov 
light (DIRC). A CsI(Tl) electromagnetic calorimeter 
(EMC) provides photon detection, electron identifica- 
tion, and 7T°, rj, and K {) L reconstruction. Finally, the 
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instrumented flux return (IFR) of the magnet allows dis- 
crimination of muons from pions and detection of K® 
mesons. For the first 214 fb -1 of data, the IFR was com- 
posed of a resistive plate chamber system. For the most 
recent 212 fb -1 of data, a portion of the resistive plate 
chamber system has been replaced by limited streamer 
tubes [20]. 



III. VERTEX RECONSTRUCTION 

In the reconstruction of the Bqp vertex, we use all 
charged daughter tracks. Daughter tracks that form a 
K® are fit to a separate vertex, with the resulting parent 
momentum and position used in the fit to the Bcp ver- 
tex. The vertex for the Stag decay is constructed from 
all tracks in the event except the daughters of Bqp- An 
additional constraint is provided by the calculated S tag 
production point and three-momentum, with its associ- 
ated error matrix. This is determined from the knowl- 
edge of the three momentum of the fully reconstructed 
Bcp candidate, its decay vertex and error matrix, and 
from the knowledge of the average position of the e + e~ 
interaction point and T(AS) average boost. In order to 
reduce bias and tails due to long-lived particles, K ° and 
A candidates are used as input to the fit in place of their 
daughters. In addition, tracks consistent with photon 
conversions (7 — ► e+e~) are excluded from the fit. To re- 
duce contributions from charm decay products that bias 
the determination of the vertex position the tracks with 
a vertex x 2 contribution greater than 6 are removed and 
the fit is repeated until no track fails the \ 2 requirement. 
We obtain At from the measured distance Az between 
the Bcp and S tag vertex with the relation Az ~ (3jcAt. 

Because there are no charged particles present at the 
S° — > n°Kg decay vertex, the ir°K^ vertex reconstruc- 
tion differs significantly from that of the u>K® and rj'K 
analyses. In ir°K® we identify the vertex of the Bcp us- 
ing the single K® trajectory from the tt + tt~ momenta and 
the knowledge of the average interaction point (IP) [21], 
which is determined several times per hour from the spa- 
tial distribution of vertices from two track events. The 
average transverse size of the IP is 180 (im x 4 [im. We 
compute At and its uncertainty with a geometric fit to 
the T(4S) — > S°S° system that takes this IP constraint 
into account. We further improve the accuracy of the 
At measurement by constraining the sum of the two B 
decay times (tcp + ttag) to be equal to 2r (r is the mean 
B° lifetime) with an uncertainty V2t, which effectively 
improves the determination of the decay position of the 
T(45). We have verified in a full detector simulation that 
this procedure provides an unbiased estimate of At. 

The estimate of the uncertainty on At for each tt°K® 
event reflects the strong dependence of the At resolution 
on the K® flight direction and on the number of SVT 
layers traversed by the decay daughters. When both 
pion tracks are reconstructed with information from at 
least the first three layers of the SVT in the coordinate 



along the collision axis (axial) as well as on the transverse 
plane (azimuthal), we obtain At with resolution compa- 
rable to that of the coK® and rj'K analyses. The average 
At resolution in these modes is about 1.0 ps. Events for 
which there is axial and azimuthal information from the 
first three layers of the SVT and for which At and the 
error on At satisfy |At| < 20 ps and (J\ t < 2.5 ps are 
classified as "good" (class g), and their At information 
is used in the time dependent part of the likelihood func- 
tion (Eq. 10). About 60% of the events fall in this class. 
Otherwise events are classified as "bad" (class b). Since 
C f can also be extracted from flavor tagging information 
alone, events of class b contribute to the measurement of 
Cf (Eq. 11) and to the signal yield in the tt°K® analysis. 

In uiKg and ij'K decays, the determination of the B 
decay vertex is dominated by the charged daughters of 
the uj and 77', so we do not require information in the first 
three SVT layers from K® daughter pions for events in 
class g. Also, since about 95% of events in these modes 
are of class g, the precision of the measurement of Cf is 
not improved by including class b events. We maintain 
simplicity of these analyses by simply rejecting class b 
events. 



IV. FLAVOR TAGGING AND At RESOLUTION 

In the measurement of time-dependent CP asymme- 
tries, it is important to determine whether at the time of 
decay of the S tag the Be p was a S° or a S°. This 'flavor 
tagging' is achieved with the analysis of the decay prod- 
ucts of the recoiling S tag meson. Most B mesons decay 
to a final state that is flavor specific; i.e., only accessible 
from either a S° or a S°, but not from both. The pur- 
pose of the flavor tagging algorithm is to determine the 
flavor of Stag with the highest possible efficiency e and 
lowest possible probability w of assigning a wrong flavor 
to i?tag- The figure of merit for the performance of the 
tagging algorithm is the effective tagging efficiency 

Q = e(l-2w) 2 , (2) 

which is approximately related to the statistical uncer- 
tainty a in the coefficients S and C through 



It is not necessary to reconstruct Stag fully to determine 
its flavor. We use a neural network based technique [22] 
to exploit signatures of B decays that determine the fla- 
vor at decay of the S tag . Primary leptons from semilep- 
tonic B decays are selected from identified electrons and 
muons as well as isolated energetic tracks. The charges 
of identified kaon candidates define a kaon tag. Soft pi- 
ons from D* + decays are selected on the basis of their 
momentum and direction with respect to the thrust axis 
of Stag- Based on the output of this algorithm, can- 
didates are divided into seven mutually exclusive cate- 
gories. These are (in order of decreasing signal purity) 
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TABLE I: Efficiencies e, average mistag fractions w, mistag 
fraction differences Aw = w(B°) — w(B°), and effective tag- 
ging efficiency Q = e(l — 2w) 2 for each tagging category from 
the Bflav data. 



Category 


e(%) 


w (%) 


Aw (%) 


Q (%) 


Lepton 


9.0 ± 0.1 


2.8 ± 0.3 


0.3 ± 0.5 


8.0 ± 0.1 


Kaon I 


10.8 ± 0.1 


5.3 ± 0.3 


-0.1 ± 0.6 


8.7 ± 0.1 


Kaon II 


17.2 ± 0.1 


14.5 ± 0.3 


0.4 ± 0.6 


8.7 ± 0.2 


Kaon-Pion 


13.7 ± 0.1 


23.3 ± 0.4 


-0.7 ± 0.7 


3.9 ± 0.1 


Pion 


14.2 ± 0.1 


32.5 ± 0.4 


5.1 ± 0.7 


1.7 ± 0.1 


Other 


9.5 ± 0.1 


41.5 ± 0.5 


3.8 ± 0.8 


0.3 ± 0.0 


All 


74.4 ± 0.1 






31.2 ± 0.3 



Lepton, Kaon I, Kaon II, Kaon-Pion, Pion, Other, and 
Untagged. 

We apply this algorithm to a sample of fully re- 
constructed, self-tagging, neutral B decays (Bfl av sam- 
ple). We use B decays to D^~(tt + , p + , a+) to mea- 
sure the tagging efficiency e, mistag rate w, and the dif- 
ference in mistag rates for B° and B° tag-side decays 
Aw = w(B°) — w(B°). The results are shown in Table I. 
The Untagged category of events contains no flavor in- 
formation and therefore carries no weight in the time- 
dependent analysis. The total effective tagging efficiency 
Q for this algorithm is measured to be (31.2 ± 0.3)%. 

Including the effects of the mistag rate, Eq. 1 becomes 

e -|A*|/T 

F(At) = — {1=fAw± (4) 

4t 

(1 - 2w) [-rifSf sin(Am d At) - C f cos(Am d At)]}. 

Finally, to account for experimental At resolution, we 
convolve Eq. 4 with a resolution function, the parame- 
ters of which we obtain from fits to the Bflav sample. 
The At resolution function is represented as a sum of 
three Gaussian distributions with different widths. For 
the core and tail Gaussians, the widths are scaled by a At- 
In addition we allow an offset for the core distribution in 
the hadronic tagging categories (Sec. IV) separate from 
that of the Lepton category, to allow for a small bias 
of At from secondary D-meson decays; a common offset 
is used for the tail component. The third Gaussian (of 
fixed 8 ps width) accounts for the few events with incor- 
rectly reconstructed vertices. Identical resolution func- 
tion parameters are used for all Bqp modes, since the 
-Stag vertex precision dominates the At resolution. 

Events without reliable At information (class b) are 
sensitive to the parameter Cf and are used to constrain 
this parameter in the n K® analysis. Integrating Eq. 4 
over At we get 

F C = \ {lT[Aw + C f (l-2w)/(l + Amy)]}. (5) 

We also account for the asymmetry in tagging efficiency 
for B° and B° decays, but, for simplicity, we assume the 
asymmetry is zero in the above equations. 



V. EVENT RECONSTRUCTION AND 
SELECTION 

We choose event selection criteria with the aid of a de- 
tailed Monte Carlo (MC) simulation of the B production 
and decay sequences, and of the detector response [23]. 
These criteria are designed to retain signal events with 
high efficiency while removing most of the background. 

We reconstruct the Bcp candidate by combining the 
four-momenta of the two daughter mesons, with a ver- 
tex constraint. The B-daughter candidates are recon- 
structed with the following decays: ir° — ► 77; rj — > 77 
(?7 77 ); 77 tt+tt-tt (?7 37r ); 7/ 77 77 7r+7r- (?^ (77)7rTr ); 

i -> Tl^ir+n- (v' n{37r)7T7r ); V' -» P°7 (V'fry), where p° -» 
7T+7T-; lu -> tt+tt-tt ; and K a s -» tt+tt- {K° +n _). In 
the 77' K® analysis we also reconstruct K® via its decay 
to two neutral pions (K° 07r0 ). The requirements on the 
invariant masses of these particle combinations are given 
in Table II. We consider as photons energy depositions in 
the EMC that are isolated from any charged tracks, carry 
a minimum energy of 30 MeV, and have the expected lat- 
eral shower shapes. 

The five final states used for B° -> t]'K° s 

rf , , K° o, and n'K° Q „. For the B° -» H 'K° 

f rj{fj)7T7r 7T U 7T U ' 'P i 7T U 7T U ' L 

channel we reconstruct the ri' in two modes: ri' , N 
and 77'/, s . Large backgrounds to the final states 
rj' pi K®, uuK®, and ivK® 07r0 preclude these modes from 
improving the precision of the measurement of CP 
parameters in rj'K and uoK°; tt°K® and ir°K® 0lT o events 
lack the minimum information for reconstruction of the 
decay vertex. 

For decays with a K® — > tt+tt - candidate we perform a 
fit of the entire decay tree which constrains the K® flight 
direction to the pion pair momentum direction and the 
production point to the Bqp vertex determined as 
described in Sec. III. In this vertex fit we also constrain 
the 77, 77', and tt° candidate masses to world-average val- 
ues [25] , since these resonances have natural widths that 
are negligible compared to the resolution. Given that the 
natural widths of the u> and p mesons are comparable to 
or greater than the detector resolution, we do not im- 
pose any constraint on the masses of these candidates; 
constraining the mass of the K° does not improve deter- 
mination of the vertex. In the tvK® and rj'K° s analyses, 
we require the \ 2 probability of this fit to be greater 
than 0.001. We also require that the K° flight length 
divided by its uncertainty be greater than 3.0 (5.0 for 
B° -» n°K° s ). 

Signal if ° candidates are reconstructed from clusters 
of energy deposited in the EMC or from hits in the IFR 
not associated with any charged track in the event [24]. 
Taking the flight direction of the K® to be the direction 
from the B° decay vertex to the cluster centroid, we de- 
termine the K® momentum Pk° l from a fit with the B° 
and K® masses constrained to world-average values [25] . 
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TABLE II: Selection requirements on the invariant masses 
of candidate resonant decays and the laboratory energies of 
photons from the decay. 

State Invariant mass (MeV) E(j) (MeV) 

Prompt n" 110 < 771(77) < 160 > 50 

Secondary tt 120 < 777(77) < 150 > 30 

?? 77 490 < 777(77) < 600 > 50 

t] 3 k 520 < m(7T- + 7r~7r") < 570 

r)'^ 945 < m(n + n~r)) < 970 

rj'pj 930 < m(TT + TY-j) < 980 > 100 

cu 735 < m(7T + 7r-7r°) < 825 

p° 470 < m{TT + Tv-) < 980 

*f° +7r - 486 < m(7r+7r-) < 510 

^"o^o 468 < m(7rV") < 528 



A. Kinematics of Y(4S) —> BB 

In this experiment the energy of the initial e + e~ state 
is known within an uncertainty of a few MeV. For a final 
state with two particles we can determine four kinematic 
variables from conservation of energy and momentum. 
These may be taken as polar and azimuthal angles of the 
line of flight of the two particles, and two energy mo- 
mentum, or mass variables, such as the masses of the two 
particles. In practice, since we fully reconstruct one B 
meson candidate, we make the assumption that it is one 
of two final-state particles of equal mass. We compute 
two largely uncorrelated variables that test consistency 
with this assumption, and with the known value [25] of 
the B-meson mass. The choice of these variables depends 
on the decay process, as we discuss below. 

In the reconstruction of B° -> uK° and B° -> rfK° s 
the kinematic variables are the energy-substituted mass 



™es = y(-s + p -p B ) 2 /^o-p| ( 6 ) 
and the energy difference 

AE = E* B - ±V~s (7) 

where (E ,p ) and (Eb,Pb) are the laboratory four- 
momenta of the T(4S I ) and the Bcp candidate, respec- 
tively, and the asterisk denotes the T(4S) rest frame. 
The resolution is 3 MeV in m ES and 20-50 MeV in AE, 
depending on the decay mode. We require 5.25 < toes < 
5.29 GeV and \ AE\ < 0.2 GeV, as distributions of these 
quantities for signal events peak at the B-meson mass in 
toes and zero in AE. 

For the B° — > rj'K® channel only the direction of the 
K h momentum is measured. For these candidates m E s is 
not determined; instead we obtain AE from a calculation 
with the _B° and K\ masses constrained to world-average 
values. Because of the mass constraint on the B°, the AE 
distribution for K® events, which peaks at zero for signal, 



is asymmetric and narrower than that of K®, events; we 
require -0.01 < AE < 0.08 GeV for r)'K°. 

For the n°K^ analysis we use the kinematic variables 
TUB and TO m ; ss . The variable tub is the invariant mass of 
the reconstructed Bcp- The variable TO m ; ss is the invari- 
ant mass of the B tag , computed from the known beam 
energy and the measured Bcp momentum with m(Bcp) 
constrained to the nominal £?-meson mass rrig DG [26]. 
For signal decays, ms and TO m i ss peak at the B° mass 
and have resolutions of ~ 47 MeV and ~ 5.4 MeV, respec- 
tively; the distribution of TOb exhibits a low-side tail due 
to leakage of energy out of the EMC. To compare the 
TO m i ss resolution with the Toes resolution a factor of two 
from the approximate relation toes ~ 
should be taken into account. The beam-energy con- 
straint in TO m j ss helps to eliminate the correlation with 
tub- We select candidates within the window 5.11 < 
J^miss < 5.31 GeV and 5.13 < ms < 5.43 GeV, which 
includes the signal peak and a sideband region for back- 
ground characterization. 



B. Background reduction 

Background events arise primarily from random com- 
binations of particles in continuum e+e~ — ► qq events 
(q = u,d,s,c). For some of the decay chains we must 
also consider cross feed from £>-mcson decays by modes 
other than the signal; we discuss these in Sees. VIA and 
VI B below. To reduce the qq backgrounds we make use of 
additional properties of the event that are consequences 
of the decay. 

For the B° — > wK® and B° — > rj'K channels we define 
the angle 9t between the thrust axis [27] of the Bcp can- 
didate in the T{AS) frame and that of the charged tracks 
and neutral calorimeter clusters in the rest of the event. 
The event is required to contain at least one charged track 
not associated with the Bqp candidate. The distribution 
of | cos#t| is sharply peaked near 1 for qq jet pairs, and 
nearly uniform for S-meson decays. The requirement is 
| cos#t| < 0.9 for all modes. 

For the r( pl decays we also define the angle Q p Aec between 
the momenta of the p° daughter tt~ and of the 77', mea- 
sured in the p° rest frame. We require |cos6^ cc | < 0.9 
to reduce the combinatorial background of p° candidates 
incorporating a soft pion that are reconstructed as decays 

with 1 cos e p dcc |~i. 

For r/'K® candidates we require that the cosine of the 
polar angle of the total missing momentum in the labo- 
ratory system be less than 0.96 to reject very forward qq 
jets. We construct the missing momentum p m ; ss as the 
difference of po and the momenta of all charged tracks 
and neutral clusters not associated with the K® candi- 
date. We project p m i ss onto Pk°, an d require the com- 
ponent perpendicular to the beam line, P^° s 3 s j_, to satisfy 

PmissJ- ~ Pk° l ± > —0.8 GeV. These values are chosen to 
minimize the uncertainty on S and C in the presence of 
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background. 

The purity of the K® candidates reconstructed in the 
EMC is further improved by a requirement on the out- 
put of a neural network (NN) that takes cluster-shape 
variables as inputs. For the NN, we use the following 
eight variables: the number of crystals in the EMC clus- 
ter; the total energy deposited in the EMC cluster; the 
second moment of the cluster energy 



^2 



(8) 



where Ei is the energy deposited in the i th crystal and r% 
its distance from the cluster ccntroid: the lateral moment 



Mlat — 



Ei 



(9) 



where E refers to the most energetic crystal and E n 
to the least energetic one; the ratio 51/59 of the en- 
ergy of the most energetic crystal (51) to the sum of 
energy of the 3x3 crystal block with 51 in its cen- 
ter (59); the ratio 59/525, where 525 is the sum of 
energy of the 5x5 crystal block with 51 in its cen- 
ter; and the absolute value of the expansion coefficients 
IZ20 1 and |^42 1 of the spatial energy distribution of the 
EMC cluster expressed as a series of Zernikc polynomials 
(C): E(x,y) = V (n,m(r,<j>), where (x,y) are 

the cartesian coordinates in the plane of the calorimeter, 
(r, <j>) are the polar coordinates of the Zernikc polynomi- 
als (0 < r < 1) and n, m are non-negative integers. The 
NN is trained on MC signal events and off-peak data 
reconstructed for the Tf' v wn K^ decay mode to return 
+1 if the event is signal- like and —1 if it is background- 
like. We check the performance of the NN on an indepen- 
dent sample of MC signal events and off-peak data recon- 
structed for the i)J, h7) „^ decay mode. Using ensembles 
of simulated experiments, as discussed in Sec. VI D, we 
find that requiring the output of the NN to be greater 
than —0.2 minimizes the average statisticial uncertainty 
on 5 and C. 

For the tt K® channel we require the x 2 probability of 
the kinematic fit to be greater than 0.001. We exclude 
events in which the absolute value of cosine of the angle 
between the beam axis and the Bcp momentum in the 
T(45) frame (cos#b) is greater than 0.9. Finally, we 
apply a cut on the event shape, selecting events with 
L2/L0 < 0.55; L{ is the i th angular moment defined as 
Li = V pj x I cos 6j\ l , where 6j is the angle with respect 
to the B thrust axis of track or neutral cluster j, pj is its 
momentum, and the sum excludes the daughters of the 
B candidate. 

The average number of candidates found per selected 
event in rj'K and u>K® is between 1.08 and 1.32, depend- 
ing on the final state. For events with multiple candidates 
we choose the one with the largest decay vertex proba- 
bility for the B meson. Furthermore, in the B° — > rf K® 
sample, if several B candidates have the same vertex 



probability, we choose the candidate with the K® infor- 
mation taken from, in order, EMC and IFR, EMC only, 
or IFR only. From the simulation we find that this algo- 
rithm selects the correct-combination candidate in about 
two thirds of the events containing multiple candidates. 

For the ir channel the average number of candi- 
dates found per selected event is 1.03. In events with 
multiple candidates we choose the one with the smallest 
value of x 2 obtained from the reconstructed mass of the 
and the tt° candidates and their respective errors. 

In the rj'K analysis, we estimate from MC the frac- 
tion of events in which we misreconstruct charged daugh- 
ters of the rf (self-crossfeed events), which dominate the 
determination of the Bcp decay vertex. We find that 
(1-4)% of events are self-crossfeed, depending on the rf 
and K° decay channel. 



VI. MAXIMUM LIKELIHOOD FIT 

The selected sample sizes are given in the first col- 
umn of Table III. We perform an unbinned maximum 
likelihood (ML) fit to these data to obtain the common 
CP-violation parameters and signal yields for each chan- 
nel. For each signal or background component j, tagging 
category c, and event class t = g (Sec. Ill), we define a 
total probability density function (PDF) for event i as 

7* Cifl = ^(AiVk, ' II 2^(4) . (10) 

k 

where T is the function F(At) defined in Eq. 4 convolved 
with the At resolution function, and if = ±1 is the B 
flavor defined by upper and lower signs in Eq.l. For 
event class t = b (for the tt°K® analysis), we define a 
total PDF for event i as 



V 1 u 



n^(4) ; 



(ii) 



where F c is the function defined in Eq. 5. The set of vari- 
ables x\, which serve to discriminate signal from back- 
ground, are described along with their PDFs Qk,j(xk) m 
Sec. VI B below. The factored form of the PDF is a good 
approximation since linear correlations are smaller than 
5%, 7% and 6% in the uK°, rf 'K° and tt°K° s analyses, re- 
spectively. The effects of these correlations are estimated 
as described in Sec. VI D. 

We write the extended likelihood function for all can- 
didates for decay mode d as 



cx p ( y., ",/,.'< ,,) 



N c , t \ 



N c , t 

n 



n jfj,t e j,cPj, c ,t 



(12) 



where rij is the yield of events of component j, fjj is 
the fraction of events of component j for each event class 
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t, ej tC is the efficiency of component j for each tagging 
category c, and N c>t is the number of events of category c 
with event class t in the sample. When combining decay 
modes we form the grand likelihood C = \\Cd- 

We fix €j iC for all components except the qq background 
to £s flav , c , which are listed in Table I. For the tt°K® 
channel we assume the same e,- iC for class g and class b 
events. For uK® and ij'K we fix fj, g = 1 because we 
accept only events of class g. 

A. Model components 

For all of the decay chains we include in the ML fit a 
component for qq combinatorial background (j — qq) , in 
addition to the one for the signal. The functional forms 
of the PDFs that describe this background are deter- 
mined from fits of one observable at a time to sidebands 
of the data in the kinematic variables that exclude signal 
events. Some of the parameters of this PDF are free in 
the final fit. Thus the combinatorial component receives 
contributions from all non-signal events in the data. 

We estimate from the simulation that charmless B 
decay modes contribute less than 2% of background to 
the input sample. These events have final states differ- 
ent from the signal, but similar kinematics, and exhibit 
broad peaks in the signal regions of some observables. 
We find that the charmless BB background component 
(j = chls) is needed for the final states uiK®, r)'K® +n _, 
and r]' K® ^ ■ We account for these with a separate com- 
ponent in the PDF. Unlike the other fit components, we 
fix the charmless BB yields using measured branching 
fractions, where available, and detection efficiencies de- 
termined from MC. For unmeasured background modes, 
we use theoretical estimates of branching fractions. 

We also consider the presence of B decays to charmed 
particles in the input sample. The charmed hadrons in 
these final states tend to be too heavy to be misrecon- 
structed as the two light bodies contained in our signals, 
and their distributions in the B kinematic variables are 
similar to those for qq. However, in the event shape vari- 
ables and At they are signal-like. We have found that 
biases in the fit results are minimized for the modes with 
i]' pj by including a component specifically for the B de- 
cays to charm states (j = chrm). Finally, for rj'r,^)™ 
we divide the signal component into two categories for 
correctly reconstructed and self-crossfeed events; we fix 
the fraction of the self-crossfeed category to the value 
obtained from MC. 



B. Probability density functions 

The set of variables Xk of Eqs. 10 and 11 is defined for 
each family of decays as 

• B a ->uK° s : {m ES , AE , F , m(ir+ ir- ir°) = m 0J ,H}, 

. B°^r,'K° s : {m ES ,A£,f}, 



. B^iKl: {AE,T}, 

• B°^>n°K°: {m B , m miss , L 2 /L a , cos 9 B }- 

Here T is a Fisher discriminant described below, and TL is 
the cosine of the polar angle of the normal to the u> decay 
plane in the to helicity frame, which is defined as the ui 
rest frame with polar axis opposite to the direction of the 
B. From Monte Carlo studies we find that including the 
rj' mass, p mass, p helicity, or w Dalitz plot coordinates 
does not improve the precision of the measurements of S 
and C. 

In Fig. 1 we show PDFs for the signal and qq com- 
ponents for the ujK® analysis, which are similar to those 
for the 7]' K® analysis. We parameterize the PDFs for 
the signal component using simulated events, while the 
background distributions are taken from sidebands of the 
data in the kinematic variables that exclude signal events. 
The parameters used in the PDFs are different for each 
mode. 

For the background PDF shapes we use the following: 
the sum of two Gaussians for <2 s ig(wEs) an( ^ Qsig(AE); 
a quadratic dependence for Q q q(AE), Q^^AE), and 
Qq-qims)] and the sum of two Gaussians for Q c hi s (A_E). 
For Qqq(mEs) and Q g g(m m i ss ) we use the function 

f(x) =xVl-x^exp[-ai-x 2 )], (13) 

with x = 2m,Es/\/s (2™ m i ss /y / s for Tr a K®) and £ a 
free parameter [28], and the same function plus a Gaus- 
sian for Qchrm("iEs) and Q c his("iEs)- For Q s i g (m B ) and 
Ssi g (w m i ss ) we use the function 

f(x) = exp ( 2 -fr-") 2 ] , (14) 

where p is the peak position of the distribution, <jl,r are 
the left and right widths, and oli,,r are the left and right 
tail parameters. For Q q q(AE) in the r)'K° analysis, we 
use the function 

f(x) =z(l-z)- 2 expo's] (15) 

where x = AE — (A.E) min , with (AE) min fixed to -0.01, 
and £' is a free parameter. 

To reduce qq background beyond that obtained with 
the cos 9t requirement described above for uK® and 
r]'K° (and the 9b and L2/L0 requirements for ir°K®), we 
use additional event topology information in the ML fit. 
The variables used include 6b, L 0} £2, and the angle with 
respect to the beam axis in the T(4S) frame of the signal 
B thrust axis (9s). For the ir°K® analysis, we use cos 9b 
and the ratio L2/L0 directly in the fit, parameterized by 
a second-order polynomial and a seven-bin histogram, 
respectively. The parameters of the L2/L0 PDF depend 
on the tagging category c in the signal component. In 
Fig. 2 we show the PDFs for signal (background) super- 
imposed on the distribution for data where background 
(signal) events are subtracted using an event weighting 
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FIG. 1: PDFs for uK° s ; from top to bottom AE, m E s, T, 
oj mass, TL, and At. In the left column we show distribu- 
tions from signal Monte Carlo; in the right column we show 
distributions for the qq component, which are taken from side- 
bands of the data in the kinematic variables that exclude sig- 
nal events. 



technique [29]. The bin widths of the L2/L0 histogram 
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FIG. 2: Distribution of (a) m m i ss , (b) ms, (c) L2/L0, (d) 
cos#s, for signal (background-subtracted) events in data 
(points) from the B° —> K^tt sample. The solid curve rep- 
resents the shape of the signal PDF, as obtained from the 
ML fit. The insets show the distribution of the data, and the 
PDF, for background (signal-subtracted) events. 



have been adjusted to be coarser where the background 
is small to reduce the number of free parameters of the 
PDF. 

For the other decay modes we construct a Fisher dis- 
criminant J 7 , which is an optimized linear combination 
of Lq, L2, |cos0 B |, and | cos 9s |- For the K\ modes we 
also use the continuous output of the flavor tagging al- 
gorithm as a variable entering the Fisher discriminant. 
The coefficients used to combine these variables are cho- 
sen to maximize the separation (difference of means di- 
vided by quadrature sum of errors) between the signal 
and continuum background distributions of J 7 , and are 
determined from studies of signal MC and off-peak data. 
We have studied the optimization of T for a variety of 
signal modes, and find that a single set of coefficients is 
nearly optimal for all. 

The PDF shape for T is an asymmetric Gaussian with 
different widths below and above the peak for signal, plus 
a broad Gaussian for qq background to account for a 
small tail in the signal T region. The background peak 
parameter is adjusted to be the same for all tagging cat- 
egories c. Because T describes the overall shape of the 
event, the distribution for BB background is similar to 
the signal distribution. 

For Qsig(^w) we use the sum of two Gaussians; for 
Qqq{ m uj) and Qchis(w w ) the sum of a Gaussian and a 
quadratic. For Q slg (Tl) and Q q q(TC) we use a quadratic 
dependence, and for Qchis(W) a fourth-order polynomial. 
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As described in Sec. IV, the resolution function in Tj (At) 
is a sum of three Gaussians for all fit components j. For 
qq background we use the same functional form Tj (At) as 
for signal, but fix the B lifetime r to zero so that Tj (At) 
is effectively just the resolution model. For the signal and 
BB background components we determine the parame- 
ters of Qk,j(x\) from simulation, and the qq background 
parameters are free in the final fit. For the signal reso- 
lution function we fix all parameters to values obtained 
from the Bfl av sample; we obtain parameter values from 
MC for the charm and charmless BB resolution models; 
we leave parameters of the At resolution model for qq 
free in the final fit. 

For the u>K® and r]'K° analyses, we use large con- 
trol samples to determine any needed adjustments to the 
signal PDF shapes that have initially been determined 
from Monte Carlo. For toes and AE in u>K® and rfK®, 
we use the decay B~ — ► n~D° with D° — > K~tt + tt , 
which has similar topology to the modes under study 
here. We select this sample by making loose require- 
ments on toes and AE, and requiring for the D° can- 
didate mass 1845 < mr> < 1885 MeV. We also place 
kinematic requirements on the D and B daughters to 
force the charmed decay to look as much like that of a 
charmless decay as possible. These selection criteria are 
applied both to the data and to MC. For T , we use a 
sample of B + — > ri'K + decays selected with require- 
ments very similar to those of our signal modes. From 
these control samples, we determine small adjustments 
(of the order of few MeV) to the mean value of the sig- 
nal AE distribution. The means and widths of the other 
distributions do not need adjustment. 

For the ui mass line shape, we use u) production in the 
data sidebands. The means and resolutions of the invari- 
ant mass distributions are compared between data and 
MC, and small adjustments are made to the PDF param- 
eterizations. These studies also provide uncertainties in 
the agreement between data and MC that are used for 
evaluation of systematic errors. For the n°Kg analysis, 
we apply no correction to the signal PDF shapes, but 
we evaluate the related systematic error as described in 
Sec. VIII. 



C. Fit variables 

For the tuK® analysis we perform a fit with 25 free 
parameters: S, C, signal yield, continuum background 
yield and fractions (6), and the background PDF param- 
eters for At, toes, AE 1 , T, to w , and H (15). For the 
five rj'K® channels we perform a single fit with 98 free 
parameters: S, C, signal yields (5), v( pi K Q s charm BB 
background yields (2), continuum background yields (5) 
and fractions (30), and the background PDF parameters 
for At, toes, AE, and T (54). Similarly the two r}'K° 
channels are fit jointly, with 34 parameters: S, C, sig- 
nal yields (2), background yields (2), fractions (12), and 
PDF parameters (16). For the tt°K® analysis we per- 



form a fit with 36 free parameters: 5*, C, signal yield, 
the means of TO m ; ss and ms signal PDFs, background 
yield, background PDF parameters for At, tub, w m i ss , 
L 2 /L , cos b (16), background tagging efficiencies (12), 
and the fraction of good events (2). For the signal, charm 
BB, and charmless BB components the parameters r 
and Ato^ are fixed to world-average values [25]; for the 
BB components S and C are fixed to zero and then varied 
to obtain the related systematic uncertainty as described 
below; for the qq model r is fixed to zero. 



D. Fit validation 

We test the fitting procedure by applying it to ensem- 
bles of simulated experiments with qq and BB charmed 
events drawn from the PDF into which we have em- 
bedded the expected number of signal and BB charm- 
less background events (with the expected values of S 
and C) randomly extracted from the fully simulated 
MC samples. We find biases (measured— expected) for 
SuK°i S, q , K a, C, q , K a, S v , k o l , and C v , K o of 0.034 ±0.010, 
0.006 ± 0.006, -0.008 ± 0.005, -0.022 ± 0.014, and 
— 0.013 ± 0.009, respectively. These small biases are due 
to neglected correlations among the observables, contam- 
ination of the signal by self-crossfeed, and the small sig- 
nal event yield in wK®. We apply additive corrections to 
the final results for these biases. For C uK o, S n o K o, and 
C v ox° we make no correction but assign a systematic 
uncertainty as described in Sec. VIII. 



VII. FIT RESULTS 

Results from the fits for the signal yields and the CP 
parameters Sf and C/ are presented in Table III. In 
Figs. 3-9, we show projections onto the kinematic vari- 
ables and At for subsets of the data for which the ratio 
of the likelihood to be signal and the sum of likelihoods 
to be signal and background (computed without the vari- 



TABLE III: Results of the fits. Signal yields quoted here 
include events with no flavor tag information. Subscripts for 
r,' decay modes denote (1) ^ (77)7r7r , ( 2 ) n' P1 , and (3) r^ (37r)7rlr . 



Mode 


# events Signal yield 


-VfSf 


Cf 




17422 


163 ± 18 


U.OO_ 29 


-0 52 +u -^ 




1470 


472 ± 24 


0.70 ±0.17 


-0.17 ±0.11 




22775 


1005 ± 40 


0.46 ±0.12 


-0.13 ±0.09 




513 


171 ± 14 


0.76 ± 0.26 


0.05 ±0.20 




1056 


105 ± 13 


0.51 ±0.34 


-0.19 ±0.30 




27057 


206 ± 28 


0.26 ±0.33 


0.04 ±0.26 




52871 


1959 ± 58 


0.53 ±0.08 


-0.11 ±0.06 




18036 


386 ± 32 


0.75 ± 0.22 


0.02 ±0.16 


VsK° L 


6213 


169 ± 21 


0.87 ±0.30 


0.19 ±0.25 




24249 


556 ± 38 


0.82l!j- 


09 +u i3 

U - UB -0.14 




21412 


556 ± 32 


0.55 ±0.20 


0.13 ±0.13 
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FIG. 3: Distributions for l}K® projected (see text) onto (a) 
mES and (b) AE. The solid lines show the fit result and the 
dashed lines show the background contributions. 



able plotted) exceeds a mode-dependent threshold that 
optimizes the statistical significance of the plotted signal. 
In ujK® the fraction of signal events with respect to the 
total after this requirement has been applied is ~ 70%, 
while in r/'K® and rfK®, the fraction of signal events is in 
the (42 — 85)% and (22 — 55)% range respectively, depend- 
ing on the decay mode. In Fig. 3 we show the projections 
onto 77t.es and AE for the wK® analysis; in Fig. 4 we show 
the projections onto tties and AE for rj'K®; in Fig. 5 we 
show the AE projections for r]'K° L . The corresponding 
information for n is conveyed by the background- 
subtracted distributions for raj and m m i ss in Fig. 2. 




m ES (GeV) AE (GeV) 



FIG. 4: Distributions for (a,b) ri' rihy)Wlr K% +lc - , (c,d) 
VpiK°+k-> ( e - f ) v'r, {37r) ^K + ^, (g,h) and 
(ij) Vp-iK°0n° projected (see text) onto (tties, AE). The solid 
lines show the fit result and the dashed lines show the back- 
ground contributions. 




AE (GeV) AE (GeV) 

FIG. 5: Distributions for (a) r/ n(iyy)WV Kl and (b) t]' n ^^Kl 
projected (see text) onto AE. The solid lines show the fit re- 
sult and the dashed lines show the background contributions. 

In Figs. 6-9, we show the At projections and the asym- 
metry (Ngo — Ng ) I (N B o + Ng a ) for each final state. In 
the u)K®, r)'K®, rfK®, and Tr°K® analyses, we measure 
the correlation between Sf and Cf in the fit to be 2.9%, 
3.0%, 1.0% and -6.2%, respectively. 



A. Crosschecks 

We perform several additional crosschecks of our analy- 
sis technique including time-dependent fits for B decays 
to the final states r/^^K* , rf^K*, and T^ (3w)mr if ± 
in which measurements of S and C are consistent with 
zero. There are only small changes in the results when we 
do any of the following: fix C — or S = 0, allow S and 
C to be different for each tagging category, remove each 
of the discriminating variables one by one, and allow the 
signal resolution model parameters to vary in the fit. 

To validate the IP-constrained vertexing technique in 
tt°K®, we examine B° — » J/ipK® decays in data where 
J/ip — > or Jfy ~ * e + e~. In these events we de- 

termine At in two ways: by fully reconstructing the B° 
decay vertex using the trajectories of charged daughters 
of the J/tp and the K° mesons (standard method), or by 
neglecting the J/ip contribution to the decay vertex and 
using the IP constraint and the K® trajectory only. This 
study shows that within statistical uncertainties of order 
2% of the error on At, the IP-constrained At measure- 
ment is unbiased with respect to the standard technique 
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FIG. 6: Data and model projections for uK® onto At for 
(a) B° and (b) B° tags. We show data as points with 
error bars and the total fit function (total signal) with a 
the solid (dotted) line. In (c) we show the raw asymmetry, 
(N B o — Ng )/{N B o + Ngo) with a solid line representing the 
fit function. 



FIG. 7: Data and model projections for rj K% onto At for 
(a) B° and (b) B° tags. Points with error bars represent 
the data; the solid (dotted) line displays the total fit function 
(total signal). In (c) we show the raw asymmetry, (N B o — 
^flo)/(^fl° +-^b°); tne solid line represents the fit function. 



and that the fit values of Sj^ K o and Cj/^^o are consis- 
tent between the two methods. 



VIII. SYSTEMATIC UNCERTAINTIES 

A number of sources of systematic uncertainties affect- 
ing the fit values of Sf and Cf have been considered. 

We vary the parameters of the signal PDFs that are 
kept fixed in the fit within their uncertainty and take 
as systematic error the resulting changes of Sf and Cf. 
These parameters include r and Am^, the mistag pa- 
rameters w and Aw, the efficiencies of each tagging cat- 
egory, the parameters of the resolution model, and the 
shift and scale factors applied to the variables related to 
the B kinematics and event shape variables that serve 
to distinguish signal from background. The deviations 
of Sj and Cf for rj'K® and r]'K® for variations of r and 
Amd are less than 0.002. 



For the Tr°K° channel as an additional systematic error 
associated with the shape of the PDFs we also use the 
largest deviation observed when the parameters of the 
individual PDFs are free in the fit. 

As a systematic uncertainty related to the fit bias on 
Sf and Cf we assign the statistical uncertainty on the 
bias obtained from simulated experiments during the fit 
validation. As explained in Sec. IV, we obtain parame- 
ters of the signal resolution model from a fit to the -Bfw 
sample instead of from a fit to signal MC. We evaluate 
the systematic uncertainty of this approach with two sets 
of simulated experiments that differ only in the values of 
resolution model parameters (one set with parameters 
from the -Bflav sample and one set with parameters from 
MC). We take the difference in the average Sf and Ct 
from these two sets of experiments as the related system- 
atic error. 

We evaluate the impact of potential biases arising from 
the interference of doubly Cabibbo-suppressed decays 
with the Cabibbo-favored decays on the tag-side of the 
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FIG. 8: Data and model projections for B° — > r\ K\ onto At 
for B° (a) and B° (b) tags. Points with error bars represent 
the data; the solid (dotted) line displays the total fit function 
(total signal). In (c) we show the raw asymmetry, (N B o — 
Ng ) /(N B o + Ng ); the solid line represents the fit function. 



FIG. 9: Data and model projections for tt°Ks onto At for 
(a) B° and (b) B° tags. Points with error bars represent the 
signal where backgroud is subtracted using an event weighting 
tcchinque [29] ; the solid line displays the signal fit function. In 
(c) we show the raw asymmetry, (N B o — Ng ) / (N B o + N B o)', 
the solid line represents the fit function. 



event [30] by taking into account realistic values of the 
ratio between the two amplitudes and the relative phases. 
For uuK® and t]'K° we estimate using MC, published mea- 
surements, and theoretical predictions that conservative 
ranges of the net values for CP parameters in the BB 
background are S = [0, 0.2] and C = ±0.1 for the charm- 
less background and S — ±0.1 and C = ±0.1 for the 
charm background. We perform a fit in which we fix the 
parameters to these values and take the difference in sig- 
nal CP parameters between this fit and the nominal fit 
as the systematic error. 

For the coK® and r/'K channels we also vary the 
amount of the charmless BB background by ±20%. For 
tt°K° we do not include a BB background component in 
the fit but we embed BB background events in the data 
sample and extract the peaking background from the ob- 
served change in the yield. We use this yield to estimate 
the change in S and C due to the CP asymmetry of the 
peaking background. We also measure the systematic er- 



ror associated with the vertex reconstruction by varying 
within uncertainties the parameters of the alignment of 
the SVT and the position and size of the beam spot. 

We quantify the effects of self-crossfeed events in the 
i]'K° analysis. For rj'K® we perform sets of simulated 
experiments in which we embed only correctly recon- 
structed signal events and compare the results to the 
nominal simulated experiments (Sec. VI D) in which we 
embed both correctly and incorrectly reconstructed sig- 
nal events. We take the difference as the systematic un- 
certainty related to self-crossfeed. For the v'^^w^^i. 
analysis, in which we include a self-crossfeed component 
in the fit, we perform a fit in which we take parameter 
values for the self-crossfeed resolution model from self- 
crossfeed MC events instead of the nominal i?fl av sample. 
We take the difference of the results from this fit and the 
nominal fit as the self-crossfeed systematic for r/'K®. The 
effects of self-crossfeed are negligible for uiK® and tt^K®. 
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TABLE IV: Summary of systematic uncertainties affecting Sf and C/. 



Source 






°ri'K° 












Variation of PDF parameters 


0.012 


0.019 


0.006 


0.009 


0.009 


0.007 


0.010 


0.012 


Bias correction 


0.010 


0.007 


0.006 


0.005 


0.014 


0.009 


0.011 


0.001 


Interference from DCSD on tag side 


0.001 


0.015 


0.001 


0.015 


0.001 


0.015 


0.001 


0.015 


BB background 


0.009 


0.010 


0.009 


0.005 






0.005 


0.001 


Signal At parameters from Bfl av 


0.002 


0.001 


0.009 


0.015 


0.004 


0.008 


0.016 


0.011 


SVT alignment 


0.011 


0.003 


0.002 


0.003 


0.004 


0.004 


0.009 


0.009 


Beam-spot position and size 


0.000 


0.000 


0.002 


0.001 


0.004 


0.003 


0.004 


0.002 


Vertexing method 














0.008 


0.016 


Self-crossfeed 






0.004 


0.001 


0.001 


0.004 






Total 


0.021 


0.028 


0.016 


0.024 


0.018 


0.021 


0.025 


0.028 



Finally, for the tt°K® analysis we examine large sam- 
ples of simulated B° -» K° s n° and B° -» J/ipK° s decays 
to quantify the differences between resolution function 
parameter values obtained from the i?fl av sample and 
those of the signal channel; we use these differences to 
evaluate uncertainties due to the use of the resolution 
function extracted from the -Bfl av sample. We also use 
the differences between resolution function parameters 
extracted from data and MC in the B a — > J/tp K®, decays 
to quantify possible problems in the reconstruction of the 
K® vertex. We take the sum in quadrature of these errors 
as the systematic error related to the vertexing method. 

The contributions of the above sources of systematic 
uncertainties to Sf and C f are summarized in Table IV. 



IX. S AND C PARAMETERS FOR B° -f rj' K 

As noted in Sec. I, the final states rfK® and rj'K® have 
opposite CP eigenvalues, and in the SM, if ASf = 0, 
then —rjfSf = sin2/3. We therefore compute the values 
of S n > K o and C n i K o from our separate measurements with 
B° — > rj'K® and B° — ► rfK®, taking — S n , K o in combina- 
tion with S v i K o, and Cy^o with C v i K o. 

To represent the results of the individual fits, we 
project the likelihood by maximizing C (Sec. VI) at a 
succession of fixed values of Sf to obtain C{—r}fSf). We 
then convolve this likelihood with a Gaussian function 
representing the independent systematic errors for each 
mode. The product of these convolved one-dimensional 
likelihood functions for the two modes, shifted in —f]fSf 
by their respective corrections (Sec. VI D), gives the joint 
likelihood for S^k - The likelihood for C„ik° is com- 
puted similarly. Since the measured correlation between 
Sf and Cf is small in our fits (Sec. VII), we extract the 
central values and total uncertainties of these quantities 
from these one-dimensional likelihood functions. Apply- 
ing the same procedure without the convolution over sys- 
tematic errors yields the statistical component of the er- 
ror. The systematic component is then extracted by sub- 
traction in quadrature from the total error. 



X. SUMMARY AND DISCUSSION 

In conclusion, we have used samples of 121 ± 13 B° — > 
uK°, 1457 ± 43 B° -> rfK°, 416 ± 29 B° -> rfK°, and 
411 ±24 B" -» 7T°X° flavor-tagged events to measure the 
time-dependent CP violation parameters 



&vK° — 


0.55±o 29 ± 0.02 


CuK% = 


-0.52±g;|g ± 0.03 


S v 'K a = 


0.57 ±0.08 ±0.02 


C n 'K° = 


-0.08 ±0.06 ±0.02 


S-k°K° = 


0.55 ±0.20 ±0.03 


C-k°K° = 


0.13 ±0.13 ±0.03 



where the first errors are statistical and the second sys- 
tematic. These results are consistent with and super- 
sede our previous measurements [14-16]; they are also 
consistent with the world average of sin2/3 measured in 
B Q -» J/ipK° s [25]. 
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